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We investigate the oscillation of a dilute atomic gas gen- 
erated by a sudden rotation of the confining trap (scissors 
mode) . This oscillation reveals the effects of superfluidity ex- 
hibited by a Bose-Einstein condensate. The scissors mode is 
investigated also in a classical gas above T c in various col- 
lisional regimes. The crucial difference with respect to the 
supcrfluid case arises from the occurrence of low frequency 
components, which are responsible for the rigid value of the 
moment of inertia. Different experimental procedures to ex- 
cite the scissors mode are discussed. 

PACS numbers: 03.75.Fi, 05.30. Jp, 32.80.Pj 



Superfluidity is one of the most spectacular conse- 
quences of Bose-Einstein condensation and has been the 
object of extensive experimental and theoretical work in 
the past, especially in connection with the physics of liq- 
uid helium Q. Indirect signatures of superfluidity in 
trapped Bose-Einstein condensed gases are given by their 
dynamic behaviour at very low temperatures Q which 
well confirms the predictions of the hydrodynamic the- 
ory of superfluids, as well as by the occurrence of spec- 
tacular interference phenomena || which point out the 
importance of coherence effects, typical of superfluids. 
However, a direct evidence of superfluidity is still miss- 
ing in these systems. 

Important manifestations of superfluidity are associ- 
ated with rotational phenomena. These include the 
strong reduction of the moment of inertia with respect 
to the classical rigid value and the occurrence of quan- 
tized vortices m. These peculiar features are the direct 
consequence of the irrotational nature of the superfluid 
flow and have been already the object of theoretical work 
also in the case of dilute trapped gases (see for example 
H and references therein) 

The purpose of this work is to focus on the oscillatory 
behaviour exhibited by the rotation of the atomic cloud 
with respect to the symmetry axis of the confining trap 
(see figure 1) and on the corresponding superfluid effects 
caused by Bose-Einstein condensation. A similar mode, 
called the scissors mode, is well known in nuclear physics 
||, where it corresponds to the out of phase rotation of 
the neutron and proton clouds, and its recent systematic 
experimental investigation has confirmed the occur- 
rence of superfluidity in an important class of deformed 
nuclei. 

In the presence of a deformed external potential the 
restoring force associated with the rotation of the cloud 



in the x-y plane is proportional to the square of the defor- 
mation parameter e of the trap (see eq.(|l]) below). The 
mass parameter is instead fixed by the moment of inertia. 
For a superfluid system this is given by the irrotational 
value and is hence proportional to e 2 . As a consequence, 
the frequency of the oscillation approaches a finite value 
when the deformation tends to zero. Viceversa, in the 
absence of superfluidity, the moment of inertia takes the 
rigid value and the scissors mode exhibits low frequency 
components. 

Let us start our investigation by considering a Bose- 
Einstein condensed gas at equilibrium in a deformed po- 
tential of the form 
m 
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with to 2 = Wq(1 + e) and u 2 = u>q(1 — e) where e gives 
the deformation of the trap in the x-y plane. For large 
enough samples one can safely use the Thomas-Fermi 
approximation for the ground state density 



n (r) = (jjl- V cxt (r))/g 



(2) 



where the strength parameter g is related to the scatter- 
ing length a by the relation g = 4:Trh 2 a/m. We consider 
a gas initially at equilibrium in the (x',y',z) frame. At 
t = one rotates abruptly the eigenaxis of the trap from 
their initial position to (x, y, z) by a small angle -9q (see 
figure 1). As a consequence of the sudden rotation, the 
system will be no longer in equilibrium and will start os- 
cillating. If the angle 0o is not too large the oscillation 
will correspond to a rotation of the cloud in the x-y plane 
(scissors mode). The equations describing this rotation 
can be easily obtained at zero temperature starting from 
the time dependent Gross-Pitaevskii equation for the or- 
der parameter which, in the Thomas-Fermi regime, takes 
the typical form of the hydrodynamic equations of super- 
fluids il: 



— — h div(nv) = 
at 

<9v / 
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Starting from (||J) one easily obtains, in the linear 
regime, the closed set of equations 



dt 



(xy) 



v + yv x ) 
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dt 



(xv y + yv x ) = ~2uj Q (xy} 



(5) 
(6) 
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involving the relevant quadrupole variable (xy) = 
J drxyn(r, t)/N. The corresponding density profile can 
be parametrized in the form 



n(r, t) = (/i - Kxt(r) - mu>Qa(t)xy) / g 



(7) 



where a is a small time dependent coefficient related to 
(xy) by 



(xy) 



2uj 2 ,a{t) 
3muj 2 u! 2 



(8) 



and (V cx t)o is the expectation value of the harmonic po- 
tential at equilibrium. The associated velocity field is ir- 
rotational and has the form v(r, t) cx V(xy). It is worth 
noticing that only if a -C e will the parametrization (Q) 
describe a rotation of the sample with the angle fixed by 
the relation a(t) = 2e9(t). If instead a is larger than e, 
then (0) corresponds to a traditional quadrupole defor- 
mation characterized by a change of the intrinsic shape 
and the connection with the geometry of the scissors is 
lost. In order to obtain a visible signal of rotational type 
in the density (0) the deformation parameter should not 
be consequently too small. Assuming a <C e (or, equiva- 
lently, <C 1) one finds, from (||), @ and the equa- 
tion 



d 2 6{t) 
dt 2 



2lu 2 6 = 



(9) 



The solution corresponding to the chosen initial condi- 
tions 0(0) = O and 0'(O) = is = o cos(w*) and 
oscillates with frequency ui = ^/2u>o- 

Let us now consider the behaviour of the system in 
the absence of Bose-Einstein condensation. The simplest 
case is the high T, classical regime where analytic solu- 
tions are available in the framework of the Boltzmann 
kinetic equations. We will use the method of the av- 
erages, already employed to discuss the damping of the 
quadrupole oscillation || . The method provides the fol- 
lowing closed set of equations 



dt 



(xy) = (xv v + yv x ) 



dt 
d 

di 



(xv y -yv x ) = 2euj%(xy) 

(xv y + yv x ) = 2((v x v y ) - uJ^(xy)) 



(10) 

(11) 
(12) 



dt 



(v x v y ) = -ujI(xv v + yv x ) - eu 2 (xv v - yv x ) 



(V X Vy) 



(13) 



where the averages are taken here in both coordinate and 
velocity space and the collisional term has been evalu- 
ated in the linear regime using a gaussian approxima- 
tion for the distribution function The relaxation 
time entering (|l^) is fixed by the relation r = 5/(47 co n) 



where 7 co u = n(0)wthc/2 is the classical collisional rate 
with n(0) the central density of the atomic cloud, v t h = 
•J 8fegT/7rm the thermal velocity and a = 87ra 2 the elas- 
tic cross-section. Notice that collisions affect only the 
equation for the variable (v x v y ). In fact the other vari- 
ables are conserved by the elastic collisions |TJJ] . It is also 
worth noticing that the angular momentum m(xv y — yv x ) 
is not a constant of motion, due to the absence of sym- 
metry in the confining potential and is coupled, through 
cq.(|ll]), to the quadrupole variable (xy). It is finally in- 
teresting to note that eqs. ( |il|Ji"3] ), with the collisional 
term set equal to zero, exactly hold also in the case of a 
non interacting Bose or Fermi gas. In this case the ef- 
fects of quantum statistics enter only through the initial 
conditions of the corresponding dynamic variables. For 
oscillations of small amplitude the density profile of the 
classical gas corresponds to the gaussian parametrization 



n(r,t) CX e -( y oxt(r)+a(t)m^a;i,)/fcT 



(14) 



and the expectation value of (xy) is given by the same 
relation (|J) previously derived for the T = Bose- 
Einstein condensed gas. The initial conditions corre- 
sponding to the sudden rotation of the gas are given 
by (xy)t=o = e&a(x 2 + y 2 )t=o, (xv y ± yv y ) t =o = and 
(v x v y )t=o = 0. The vanishing of (xv y ± yv y )t=o follows 
from the absence of currents at t = 0, while the vanish- 
ing of (v x v v )t=o is the consequence of the initial isotropy 
of the velocity distribution. - By using the relationship 
a(t) = —2e6(t) holding for small rotational angles the 
equations of motion (JlC -11) can be usefully rewritten in 
the form 



d A 6 
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,d0\ 
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= 
(15) 



and the initial conditions take the form 0(0) = 0o, 
0"(O) = -2w 2 0(O) and 0'(O) = 0"'(O) = 0. Notice that 
eq.(|l^) differs from the corresponding equation ([)]) for 
the superfluid regime. It gives rise to different solutions 
propagating at high and low frequency. For small values 
of the deformation parameter e, the former can be iden- 
tified with the £ z = 2 irrotational quadrupole oscillation, 
i.e. the classical counterpart of the superfluid oscillation 
described by (^|). The latter instead corresponds to the 
rotational mode of the system and is absent in the super- 
fluid case. Actually, if e becomes too small, the variable 
looses its geometrical meaning. In this case the physical 
variables are on the one hand the moment (xy), which 
characterizes the high lying quadrupole mode and is cou- 
pled with the variables (xv y +yv y ) and (v x v y ), and on the 
other hand the angular momentum m(xv y — yv x ) which 
becomes a constant of motion (see eq.(ll)). 

Let us discuss the different collisional regimes pre- 
dicted by JlS|)- In the collisionless regime (first term in 



2 



(|l5|)) the two frequencies are undamped and given, re- 
spectively, by | oj x ± oj y | . In the opposite, hydrodynamic 
limit (second term) only the high-lying oscillation sur- 
vives with frequency u> = V^o, while the low- lying so- 
lution becomes overdamped as r — > 0, in agreement with 
the diffusive nature of the transverse waves predicted by 
classical hydrodynamics jOJ. A maximum damping is 
obtained for luqt ~ 1 a condition easily achievable in 
current experiments. 

The chosen initial condition, corresponding to a sud- 
den rotation of the sample, gives rise to the excitation 
of both the low and high frequency modes. The result- 
ing time evolution of the observable 9(t) is shown in fig. 2 
and 3 for two different collisional regimes and the choice 
e = 0.3. In fig. 2 we have made the "collisionless" choice 
€luqt ~ 6, so that both the high lying and low lying modes 
have small damping. The figure clearly shows the com- 
bined signals propagating with frequencies | u> x ± u> y | 
respectively. In fig. 3 we have instead made the choice 
ooqt ~ 0.15 corresponding to large damping. With such 
a choice the low frequency rotational mode is overdamped 
and the remaining oscillation exhibited by the curve is 
due to the high frequency component The achievement 
of the overdamped, "hydrodynamic" regime for the low 
frequency mode, is favoured by the choice of small values 
of e. 

From the comparison between the curves of fig. 2 and 
3 it emerges very clearly that the main feature of the 
superfluid regime (dashed line in fig. 2 and fig. 3) is the 
absence of low frequency components. This behaviour 
should be easily verificable experimentally. 

Untill now, we have discussed the superfluid case (T = 
0) and the classical regime T > T c . Below T c the sys- 
tem can be described in terms of a two-fluid model and 
one needs two different angles to describe the motion of 
the system. Neglecting interaction effects between the 
two fluids the superfluid component would be governed 
by (||) while the thermal part would evolve according to 
( |l5| ) . Of course in this case the relaxation time r should 
be evaluated by taking into account Bose statistics. In- 
clusion of interaction effects between the condensate and 
the thermal component would lead to a damping of the 
superfluid oscillation, similarly to what happens in the 
case of the quadrupole mode fll^l ■ 

The sudden rotation of the trap is not the only way to 
excite the scissors mode. In the last part of this letter 
we discuss an alternative procedure which further em- 
phasizes the superfluid nature of the Bosc-Einstein con- 
densed gas. We consider a gas initially in equlibrium 
within a trap rotating with frequency Q. Experimen- 
tal techniques to achieve rotating configurations of this 
type are in progress Jl^] . At the time t = we suddenly 
stop the rotation of the trap and the gas, due to its in- 
ertia, will start rotating. In the absence of superfluidity 
the moment of inertia is large and the oscillations will 
be characterized by low frequency components and large 



amplitudes. Viceversa, if the system is superfluid the 
moment of inertia is small (proportional to e 2 ) and the 
oscillations will be characterized by large frequencies and 
small amplitudes. The equations of motion describing the 
rotation are still given by (g) and @ but the initial con- 
ditions arc different, corresponding to the presence of a 
current term in the system at t = 0. One finds 0(0) = 
and #'(0) = Cl. In the classical case one has the addi- 
tional conditions 0"(O) = 0"'(O) = 0. The amplitude of 
the oscillations, in the collisionless regime, scale as fi/wo 
for the superfluid and as Cl/(oJoe) for the classical gas. 
The difference is due to the fact that at t = the veloc- 
ity field generated by the rotation is very different in the 
two cases. In the superfluid it is given by the irrotational 
form v = — fleV(xy), while in the classical case by the 
rotational form v = flxr. If the amplitude of the oscilla- 
tion becomes too large the dynamic variable (xy) cannot 
be longer simply connected with the angle 0. This easily 
happens in the classical case where the amplitude of the 
oscillation is magnified by the factor 1/e. 

The excitation of the scissors mode based on a sud- 
den switching off of the rotation of the trap would al- 
low also for an explicit determination of the moment 
of inertia O of the gas defined by the linear relation 
m(xv y — yv x )t=o = QQ/N. In fact the knowledge of 
the time evolution of (xy) permits to calculate directly 
the value of the angular momentum at t = 0. From 
eq.(|Tl]), which holds for classical as well as for quantum 
systems, one has m{xv y — yv x ) t= o = — 2ecjQ J du>F(u))/uj 
with F(tu) defined by 



du}F(u})(e i[ 



f )/2* 



(16) 



On the other hand, from the model independent equa- 
tion (10), one easily finds the result (xv y + yv x )t=o — 
J cLuF(uj)uj so that, using the relation (xv y + yv x )t=o = 
£l(x 2 — y 2 )t=o predicted by linear response theory pH , 
one finally obtains the useful expression 



6 = e rig (w a 



x) (x 2 + y 2 )t=o Jdu;F(uj)u 1 ' 

relating the moment of inertia of the system to the mea- 
surable Fourier signal F(u>). In this equation O r i g = 
Nm(x 2 + y 2 ) is the rigid value of the moment of inertia. 
Notice that eq. (|l7]) provides an exact relationship for the 
moment of inertia holding for classical as well as Bose or 
Fermi interacting systems confined by a harmonic trap 
of the type ([!]). If the system oscillates with the single 
frequency ui — \/2<^o> as happens in the T = Thomas- 
Fermi superfluid regime, one immediately finds the result 
O = e 2 O r ; g corresponding to the irrotational value of the 
moment of inertia. Viceversa, in a classical gas or in a 
normal Fermi gas trapped by a harmonic potential, the 
response is dominated by low frequency components, of 
order etoo, and eq.(p^) yields the rigid value for the mo- 
ment of inertia. 
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FIG. 1. Exciting the scissor mode : initially the gas is ther- 
malized in an anisotropic trap. One then abruptly rotates the 
eigenaxis of the trap by a small angle. 
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FIG. 2. Angle 8 as a function of time for a classical gas 
(full line) in the "collisionless" regime (euor — 6 and e = 0.3) 
and for a superfluid (dashed line). 
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FIG. 3. Angle 8 as a function of time for a classical gas 
(full line) in the "hydrodynamic" regime (euoT ~ 0.15 and 
e = 0.3) and for a superfluid (dashed line). 
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